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Abstract. The Virial Thoerem (VT) is a mathematical expression 
obtained from the equation of motion for a fluid, which describes the 
energy budget of particular regions within the flow. This course re¬ 
views the basic theory leading to the VT, discusses its applicability 
and limitations, and then summarizes observational results concern¬ 
ing the physical and statistical properties of interstellar clouds which 
are normally understood in terms of the VT, in particular the so- 
called “Larson’s Relations”. In all cases, the standard notions as well 
as relevant counterpoint are presented. In particular, the difficulties 
arising when the medium is fully turbulent are discussed. 
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1. Introduction 

The interstellar medium is an extremely complex mixture of gas, dust, and 
cosmic rays, threaded by a ubiquitous magnetic field, and stirred by a 
wide variety of energy sources such as supernovae, expanding HII regions 
and bipolar outflows. Progressively higher-density regions in this medium 
constitute what we normally call intercloud medium, cloud complexes, dif¬ 
fuse clouds, molecular clouds, clumps and cores, respectively, although the 
boundaries between these classifications are fuzzy (a natural consequence of 
drawing arbitrary classification boundaries in a continuum). Table gives 
the ranges of densities and temperatures for the above structures (adapted 
from Jura 1987). 

The gaseous component of this medium should be reasonably well de¬ 
scribed as a compressible, self-gravitating fluid (e.g, Shu 1992, Ch. 1) obey¬ 
ing the magnetohydrodynamic (MHD) equations, since even very low ion¬ 
ization fractions allow coupling of the motions to the magnetic field (Mestel 
&: Spitzer 1956; see also Shu 1992, Ch. 27). Prom one of these equations, 
a very important result called the Virial Theorem (VT) can be obtained, 
which describes the energy balance (or “budget”) of particular regions within 
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Table 1. Structures in the interstellar medium. 
(Adapted from Jura 1987.) 


Name 

n(cm 

r(K) 

Filling factor 

Hot intercloud 

3 X 10-3 

10® 

0.1 -0.7? 

gas 

Warm gas 

0.1 

103 - 10^ 

0.4 

Diffuse clouds 

< 10 

50 - 100 


Cirrus clouds 

10 - 103 

10 - 100 


Dark clouds 

> 103 

10 

0.01? 

GMCs 

> 103 

15-40 


H II regions 

> 10 

10^ 


SNR 

> 1 

10^ - 10^ 
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the medium. In this course, we first review the derivation of the VT (§ 

i), stressing in particular the physical meaning of each term, including the 
often-neglected surface terms, and the validity of frequently used simplifying 
assumptions (e.g., Shu 1992). In § |3] we then discuss various applications 
of the Virial Theorem, both to idealized cases and to real interstellar cloud 
properties, pointing out existing incompletenesses. Finally, in § ^ contains 
the conclusions. Throughout this course, we shall ignore the flux-loss effects 
of ambipolar diffusion (e.g, Shu 1992). 


2. The MHD Equations and the Virial Theorem 


2.1. The MHD equations 

In order for the gaseous component of the ISM to be adequately described 
by a continuum approximation, it is necessary that the mean free paths of 
the gas molecules be much smaller than the characteristic length scale of the 
systems under consideration (e.g., Currie 1974). This is generally satisfied 
in most kinds of interstellar structures, since the characteristic length scales 
are normally at least some 10® times larger than the mean free paths (e.g., 
Shu 1992). The use of the MHD equations is thus amply justified. These 
equations, describing the evolution of a magnetized, compressible gas with¬ 
out dissipative terms, are (e.g.. Cowling 1976; Spitzer 1978; Shu 1992; Shore 
1992): 

i) The Continuity Equation, expressing mass conservation: 

+ V • ipu) = 0; (1) 

ii) The Momentum Balance Equation 

- h pu ■ Vu = —VP — pVif + —(V X B) X B — 2pQ, x u; (2) 

ot dvr 

iii) The Internal Energy Balanee Equation 

dc 

-— -I- u • Ve = -(7 - l)eV • u -I- T - pA; (3) 

at 

iv) The Magnetie Eield Equation (flux freezing) 


dB 


V X (u X B); 


( 4 ) 


v) Poisson’s Equation 


= InGp. 

Additionally, we assume an ideal-gas equation of state P = (7 — l)pe. 


( 5 ) 
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In these equations, p is the mass density of the gas, u is the fluid velocity, 
P is the thermal pressure, e is the specihc internal energy, B is the magnetic 
field, and (p is the gravitational potential. Additionally, V and A symbolically 
represent any heating and cooling sources present in the medium, and are 
normally functions of the density and temperature in the MHD description 
(neglecting radiative transfer and chemical composition). 

Let us now briefly discuss the meaning of the main terms in these equa¬ 
tions. In eq. (|^, which is Newton’s Second Law in per-unit-volume form, 
the second term on the left hand side (LHS) is the so-called non-linear or 
advective term, and expresses the momentum transport by the velocity field 
itself. This is the term that describes the effects of turbulence in the medium. 
The first term in the right hand side (RHS) gives the force (per unit vol¬ 
ume) exerted by the thermal pressure gradient. The second term is the 
self-gravitational force, the third term is the Lorentz force, and the last 
term is the Coriolis force, due to Galactic rotation with an angular velocity 
Q. 

In eq. (^, the second term on the LHS expresses internal energy trans¬ 
port by the velocity field (convection). The first term on the RHS is the PdV 
work. In eq. (^), the RHS is the flux freezing condition, implying that mag¬ 
netic field lines are dragged along with the matter. Finally, eq. (|^ expresses 
that the source of the gravitational potential is the mass density. 

Important physics left out from this description are, as mentioned above, 
the radiative transfer and the chemical and ionization fraction evolution. 
These can be crudely incorporated in the equations above by means of model 
terms, such as T and A in eq. (^, or by making the equation of state vary 
in different regimes. This allows an approximate description of the effects of 
these processes on the dynamics, although it cannot allow any predictions 
on those processes themselves. Finally, note that the MHD equations hold 
in the so-called “MHD approximation” which assumes an conductivity of 
the medium. 


2.2. The Virial Theorem 

The scalar VT is obtained by dotting the momentum equation (eq. ^) with 
the position vector x and integrating over some volume V. Traditionally, the 
volume V is taken in a Lagrangian frame of reference moving with the fluid 
(e.g., Spitzer 1978; Shu 1992). An Eulerian (i.e, in a fixed frame of reference) 
version of the VT has recently been derived by McKee & Zweibel (1992, 
hereafter MZ92). This may be a more convenient form for computation in 
practice (e.g., MZ92; Ballesteros-Paredes & Vazquez-Semadeni 1997), but 
here we discuss the Lagrangian version for consistency with most treatments. 
Neglecting the Coriolis force, which is negligible at molecular cloud scales. 
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one obtains 


f ^-VPdV-j p^-VipdV+^ f x-(VxB)xBdy, (6) 
V Jv 47r Jv 


where d/dt = d/dt + u ■ V is the total or Lagrangian derivative operator. 
The LHS can be rewritten as 



where we have defined dm = pdV, and have used Reynolds’ Transport The¬ 


orem 



IV 


da 

-'m 


+ V • (au) 


dV 


( 8 ) 


to show that ^ / adm = / ^dm. Also, the second equality in eq. (@) defines 
the moment of inertia I. 

For the RHS, consider the following. First, according to Newton’s Law 
of Gravitation (or from the solution of Poisson’s equation) the gravitational 
term can be written as 


W = - p{x)x-VipdV = -G 


IV 


V Jail space 


d^i£h>h^z2hav'dv. (9) 


X — X' 


Noting that the integrand is nearly symmetric with respect to the primed 
and unprimed variables, we can replace x by (x — x')/2 and write 


W = --G 


1 r r 

2 Jv Jail space ^ | 


( 10 ) 


Thus, W is the gravitational energy of the mass distribution within volume 
V, given the mass distribution of the whole space.Q Note that in Shu (1992), 
both integrals are evaluated over volume V. This approximation is valid 
only if the mass outside volume V can be neglected (see Spitzer 1978), a 
fact which is often overlooked. 

Secondly, the remaining thermal pressure and magnetic terms can be 
dealt with using Gauss’ divergence theorem, and the facts that V • x = 3 
and that (in tensor notation) dxijdxj = 6ij. Additionally, for the magnetic 
term, it is most convenient to work with Maxwell’s magnetic stress tensor 


^Interestingly, this is not true in two dimensions (2D) (Ballesteros-Paredes & Vazquez- 
Semadeni, in preparation), due to the different spatial variation of the gravitational po¬ 
tential, analogously to the well-known situation for electromagnetic fields. This is un¬ 
fortunate, since many numerical simulations are performed in 2D due to computational 
limitations. 
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T = Tij = BiBj/An — |Bp5jj/87r, such that the flux freezing term satisfies 
(V X B) X B)/47r = V • Tij. Thus, the VT finally reads 


^ = J pu^dV +(3j PdV - f dS 


t 

2 dP 
+ 


[ f fpri^dvdv 

Jv J aW space |x X | 

= 2ii: + ( 2 c/ - ^ Px • dS) + (m + ^ X • T • ds) + IT, 


-1g 

2 


( 11 ) 


where S is the surface enclosing V. 

Several comments and warnings are in order. First, note that the LHS is 
the second time derivative of the moment of inertia of the cloud. The cloud 
is said to be in virial equilibrium if the RHS is zero. However, note that in 
particular this includes the case of a constant rate of change in the cloud’s 
moment of inertia. Thus, virial equilibrium does not necessarily imply a 
stationary or static cloud configuration, although the expression is almost 
always interpreted to mean so. 

Second, note that both the thermal pressure and the magnetic terms 
contain a “volumetric” and a “surface” contribution. In fact, if either the 
thermal pressure or the magnetic field are uniform, then the corresponding 
surface and volumetric terms cancel each other, indicating balance between 
internal and external stresses. Conversely, the surface terms can only be 
ignored if the corresponding fields can be neglected at the boundary surface. 
However, it is frequently encountered in the literature that virial balance is 
considered only among the volumetric terms. 

Third, the thermal pressure surface term describes the “confining” ef¬ 
fect of the external pressure on the cloud (volume V). A similar comment 
applies to the magnetic surface term, although in this case this term is not 
exclusively confining, as it includes the effects of both the external mag¬ 
netic pressure as well as the magnetic tension along field lines. Moreover, in 
the Eulerian version of the theorem, an analogous surface term of the form 
/ X • puu • dS for the macroscopic velocity field also appears (MZ92), where 
puu is a momentum transfer tensor reminiscent of the Reynolds’ stress ten¬ 
sor (see, e.g., Landau & Lifshitz 1987, § 42; Lesieur 1990, Ch. VI). However, 
due to the strongly fluctuating nature of the turbulent velocity field, which 
may contain large-scale motions (comparable to the cloud’s size), this term 
can hardly be considered as producing a confining effect; instead, it produces 
transport of momentum across the cloud’s boundary and a redistribution of 
the mass within the cloud, thus contributing to the change in its moment of 
inertia. Neglect of this type of exchange across an Eulerian boundary, or of 
the high mobility of a Lagrangian boundary, led MZ92 to the conclusion that 
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the intercloud velocity component perpendicular to the surface of the cloud 
must be small at small distances from the boundary surface S, a result which 
appears questionable in the light of numerical simulations of turbulence in 
the ISM (Passot et al. 1988; Leorat et al 1990; Vazquez-Semadeni et ah 
1995; Passot et al. 1995), which suggest an extremely dynamical scenario. 

The last point leads to a practical warning. The choice of the vol¬ 
ume V and its boundary S in the actual ISM or in numerically simulated 
flows constitutes a difficult task as soon as high enough resolution is avail¬ 
able. The VT has traditionally been applied to idealized spheroidal, static 
clouds either in magnetohydrodynamic equilibrium (possibly assisted by ex¬ 
ternal pressure confinement; e.g., Strittmatter 1966; Mouschovias &: Spitzer 
1976; see also Shu 1992), equilibrium between internal (micro)turbulence and 
self-gravity (Chandrasekhar 1951; Bonazzola et al. 1987; Vazquez-Semadeni 
&: Gazol 1995), or combinations thereof (e.g, Myers & Goodman 1988a,b; 
Mouschovias & Psaltis 1995). In these cases, the clouds have well defined 
boundaries. However, if the clouds are highly dynamical, even possibly with 
fractal structure (Scalo 1990; Falgarone et al. 1991), and they are observed 
with high enough resolutions (e.g.. Bally et al. 1987; see also the numerical 
results of Ballesteros-Paredes & Vazquez-Semadeni 1997), then important 
difficulties arise. In a Lagrangian description, the cloud boundaries (the sur¬ 
faces S) move in an extremely complex way, getting stretched and distorted. 
In an Eulerian description, a surface S dehned at a particular time will 
have an extremely amorphous shape, possibly being very elongated or hav¬ 
ing “tentacle”-like protrusions. Thus, at subsequent times, mass will have 
entered or left the contained volume V, rendering the identity of the cloud 
dubious (Ballesteros-Paredes & Vazquez-Semadeni 1997). In this case, it is 
likely that all the terms in the RHS of the VT have comparable importance, 
particularly the surface terms. With this caveat in mind, we now proceed to 
review some of the best known applications of the VT to interstellar clouds. 


3. Applications 
3.1. Ideal cases 

Pressure confinement. Taking B = u = 0 and neglecting self-gravity in 
eq. one gets the condition of virial equilibrium, d?I/dt^ = 0 2U = 
• dS, indicating balance between internal and external pressure. This 
is a good example of a case in which volume V can be chosen in a completely 
arbitrary way, since this is the case, for example, of any parcel of gas in the 
air in the room you are in. 

Hydrostatic equilibrium and turbulent support. Assuming B = u = 0 and 
a vanishing thermal pressure outside volume V, but keeping self-gravity, we 
obtain d?I/dt^ = 0 ^ W = —2U, indicating hydrostatic equilibrium. In- 
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eluding a turbulent velocity field, the condition of virial equilibrium becomes 
W = —2{U + K), indicating that turbulence can provide, in principle, ad¬ 
ditional support against gravity. However, note that this assumes that the 
turbulence is microscopic [microturbulence)^ since it is implicitly required 
that the velocity field itself does not alter the mass distribution within the 
volume. 


Magnetic support. Assume now a spherical cloud (a density peak) of ra¬ 
dius R, with a uniform field B = inside and a vanishing field out¬ 

side, and u = P = 0. It is easy to show that for such a cloud, W = 
—G /[M(r)/r]dM(r) = —3GM^/5R and M = B^R^/Q = cj>^/[Gir'^R), where 
M[r) is the mass interior to radius r and cf) is the magnetic flux, which, 
according to the flux-freezing condition, eq. i) , is conserved in the absence 
of magnetic dissipation. Since (j) is constant, both W and M. scale as R~^ 
upon contraction, and the ratio Wj Mi is constant. This implies that there is 
a critical mass-to-flux ratio [Ml<fl)cnt = (5/18)^/^(7r^G)“^/^ below which a 
cloud cannot collapse gravitationally, since the magnetic energy Mi is larger 
than the gravitational energy W (Strittmatter 1966). Modifications to this 
critical ratio due to the geometry of the cloud were studied analytically also 
by Strittmatter (1966) and numerically by Mouschovias &: Spitzer (1976), in 
the presence of an external pressure. Bertoldi & McKee (1992) have studied 
the stability of spheroidal, pressure-confined clumps, giving estimates of the 
Jeans and magnetic critical masses. 

Finally, let us note that throughout this course, we discuss the scalar 
version of the VT. A tensor form can also be derived, which allows the treat¬ 
ment of multidirectional transport phenomena. Zweibel (1990) has given a 
detailed account of magnetic support and stability of spheroidal, quiescent 
cloud configurations using the tensor VT. 


3.2. Properties of molecular clouds. Larson’s relations 

Molecular clouds exhibit a number of correlations between their various phys¬ 
ical properties, such as their size, velocity dispersion, magnetic field strength 
and (possibly) density and mass. Larson (1981) first noticed, from a collec¬ 
tion of data from several surveys available to him, that the velocity dispersion 
a and number density n of molecular clouds scaled with their size as 


^km s 


1.1 


/ R \0.38 
Vpc7 ’ 



( 12 ) 


A large number of subsequent studies (e.g.. Torrefies et al. 1983; Dame et al. 
1986; Falgarone &: Perault 1987; Myers &: Goodman 1988a; Falgarone et al. 
1992; Miesch & Bally 1994; Wood et al. 1994; Caselli & Myers 1995) have 
found similar scalings, although the currently most favored forms are 


a ~ 


(13) 
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and 

p^R-^ (14) 

where p = mn is the mass density, and m is the mean molecular mass. These 
scaling relations have been interpreted, since their discovery, as evidence 
for virial equilibrium of molecular clouds. Larson (1981) noted that, as a 
consequence of these relations, the ratio 2GM/a‘^R was roughly constant for 
his cloud sample, indicating virial balance between the turbulent velocity 
dispersion and self-gravity {W ~ K). A problem with this result however, 
was that the observed linewidths imply in general highly supersonic velocity 
dispersions. Therefore, it was expected that such random motions should 
produce shocks and dissipate rapidly, becoming subsonic (e.g., Mestel 1965a, 
b). This problem disappears, however, if one introduces magnetic fields into 
the picture, since Alfven waves do not dissipate as rapidly, and additionally 
for clouds with nearly the critical mass-to-flux ratio, the Alfven velocity 
nearly equals the virial velocity (Shu et al. 1987; Mouschovias 1987). That 
is, from IT ~ A4 one gets 


4:Trp 


= vl 


rsj 


GM 


~ GpR^, 


(15) 


while from W ~ iL, we obtain 


2 

(j ~ 


GpR^ 


(16) 


where ua is the Alfven speed^ This in principle could solve the dissipation 
problem, since in this case velocity d ispe rsions a up to the Alfven speed 
do not induce shocks (although see 


while explaining the observed 
virial value of the velocity dispersion. If additionally one assumes that for 
an ensemble of clouds the density-size relation p ~ R~^ is satisfied, then 
(t2 ~ GpR^/R r\j R, automatically satisfying the velocity dispersion-size 
relation. That is, out of the two Larson’s relations and the virial equilibrium 
condition, only two are independent. Furthermore, the relation p ~ R~^ can 
be interpreted as a consequence of the clouds being near the critical mass- 
to-flux ratio at roughly constant magnetic held strengths. Indeed, from eq. 
(15) it follows that p ~ R~^ for constant B. 

Myers & Goodman (1988a) tested these ideas on a cloud sample with 
existing magnetic held strength determinations from Zeeman measurements. 
Assuming cr ~ ua, it follows from eqs. (|T^ and (0) that the held strength 
necessary for virial equilibrium is i?eq ~ jR. Myers & Goodman (1988a) 
plotted the observed magnetic held strength versus i?eq; hnding the two were 
equal to within factors of 3 (hg. I). 


■^Xie (1997) has pointed out that it is actually the magnetic fluctuations AS that should be 
considered in these expressions, rather than the mean fleld B. However, the two generally 








Figure 1. Observed vs. virial equilibrium values of the magnetic 
field for the cloud sample of Myers & Goodman 1988a. 






Figure 2. Observed magnetic field strength vs. virial density 
n = 15(327rln2mG)“^(Au/i?)^ for the cloud sample of Myers & 
Goodman 1988a. The solid lines indicate the virial equilibrium 
model for the indicated velocity dispersions. 

Additionally, they also found that i?obs ~ in agreement with ~ 

/p (again derived from combining relations and j^), since in their 

sample a had relatively little scatter compared with the other quantities. 
The resulting correlation had a scatter consistent with the scatter in a (fig. 
§). Note, however, that in the latter plot, the “virial” density they used 
was an estimate assuming virial equilibrium {p oc (Au/i?)^), not an actual 
observed value, and thus this plot can only be taken as a proof of self- 
consistency of the virial equilibrium hypothesis. 

Myers &: Goodman (1988b) also considered the fact that the observed 
linewidths contain both a thermal and a nonthermal component, explaining 
deviations from Larson’s (1981) trends at very small scales (cloud cores of 
sizes ~ 0.1 pc). According to relation (|I^), which refers to the macroscopic, 
nonthermal velocity dispersion, at small enough cloud sizes the nonthermal 
component becomes smaller than the thermal (the velocity dispersion be¬ 
comes subsonic). Indeed, Myers & Goodman (1988b) find a drop in the 
ratio of nonthermal to gravitational energy in their cloud sample 

for small cloud sizes. It is important to note, however, that the trend is 
marginal when the raw data are used, but clearly noticeable only when the 
assumption of virial equilibrium is introduced (replacing the true gravita¬ 
tional energy obtained from direct measurements of the size and density by 
the total kinetic (thermal -|- nonthermal) energy). From this result, Myers 


turn out to be of comparable magnitude (Myers & Goodman 1991), and use of B is 
justified, alhough the difference should be kept In mind. 
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& Goodman (1988b) suggest that the uncertainties in arise more 

from uncertainties in the density and size than from uncertainties in velocity 
dispersion and tempertaure. However, an obvious alternative possibility is 
that the scatter is real, indicating departures from precise virial equilibrium. 
Their plot “assuming virial equilibrium” is in reality just a plot of the ratio 
of nonthermal to total kinetic energy, thus necessarily being very tight at 
large scales if the velocity dispersion does increase with size while the tem¬ 
perature remains roughly constant, but it is not a test of virial equilibrium 
in the clouds. On the other hand, further supporting evidence for virial 
equilibrium has been provided by Goodman et al. (1993), who find good 
agreement between the observed and the “virial” values of /3, the ratio of 
rotational to gravitational energy in a sample of dense cores. 

The role of thermal plus nonthermal (TNT) motions has been further 
investigated by Fuller & Myers (1992) and by Caselli & Myers (1995), finding 
respectively that the scaling of the nonthermal part of the velocity differs 
from that of the total velocity dispersion, and that the scaling is different for 
massive cores than for low-mass cores. Furthermore, Fuller Sz Myers (1992) 
have distinguished between cloud-to-cloud scaling, i.e., that obtained when 
the sample consists of various clouds observed with the same tracer (sensitive 
to one particular density range) and line-to-line scaling, i.e., the scaling 
obtained when using various tracer molecules (sensitive to a variety of density 
ranges) on a single cloud. The latter can be used to test for virialization 
and the density structure of individual clouds as a function of the distance 
from the core. An even finer subdivision of possibilities has been proposed 
by Goodman et al. (1997). These authors have additionally found that 
the slope of the a-R relation using single-line, single-cloud measurements, 
appears to decrease to nearly zero at scales ~ 0.1 pc, and that the column 
density N changes its size dependence from N ~ R~^ ‘^ at scales 0.1 pc to 
N ~ R~^'^ at scales 0.1 pc. They interpret these results as evidence for 
a transition from non-thermal to predominantly thermal support, and the 
onset of “velocity coherence” at those scales. 

Finally, it should be noted that a number of other possible mechanisms, 
capable of originating Larson’s relations or either one thereof, have been 
proposed. For example, Larson (1981) himself suggested that the density- 
size relation might be due to planar shocks which keep the column density 
constant along the direction of shock propagation. Ghieze (1987) has pro¬ 
posed that an ensemble of clouds on the verge of gravitational instability 
in an environment providing external pressure should also satisfy Larson’s 
relations. In the chapter on Turbulence in Molecular Glouds, a discussion on 
mechanisms related to hydro- or magnetohydrodynamic turbulence is given. 
Other scenarios are unfortunately out of the scope of this course. 




Figure 3. Nonthermal linewidth vs. size for 14 objects from 
Myers &: Goodman 1988a. No clear correlation is seen. (From 
Mousctiovias & Psaltis 1995.) 


3.3. Counterpoint 

As it may have started to be apparent to the reader from the discussion of the 
previous section, the issue of cloud virialization as the origin of the Larson 
(1981) relations is not completely settled. Besides the caveats inherent to the 
work of Myers &: Goodman (1988a, b), which are the standard references for 
observational determinations of virial equilibrium in clouds, other somewhat 
contradictory evidence exists. 

First, let us distinguish between a cloud being in virial equilibrium and 
it satisfying the VT. Actually, all clouds, and in general all fluid parcels 
in the ISM, satisfy the VT, as it is a direct consequence of the momen¬ 
tum conservation equation. Instead, virial equilibrium means d?I/dt^ = 0 
specifically. 

Goncerning the velocity dispersion-size relation, Scalo (1990) and Mous- 
chovias & Psaltis (1995) have pointed out that the original compiled by 
Myers & Goodman (1988a) show essentially just scatter for this relation (fig. 
p.3j ). Mouschovias & Psaltis (1995) “solve the problem” by plotting instead 
the ratio vs. the magnetic field strength B. In reality, however, this is 

nothing but the same Bobs vs. Bgq plot shown by Myers &: Goodman (1988a). 
Moreover, Mouschovias &: Psaltis (1995) use this result as an argument in 
favor of the a-R relation being a consequence of magnetic support. However, 
they arrive at this conclusion by assuming B does not “vary much from 
place to place in the ISM under conditions suitable for the formation of 
self-gravitating clouds”; yet, the magnetic field strength data they consider 
vary by three orders of magnitude! In summary, the cloud data discussed 
in Myers & Goodman (1988a) are consistent with magnetic support, but do 
not satisfy the dispersion-size relation. 





14 


Enrique Vazquez-Semadeni 


Furthermore, there exists both observational and numerical evidence 
that the density-size relation is not always satisfied, implying that the col¬ 
umn density is not constant. It was already pointed out by Larson (1981) 
himself that the density-size relation might not be a real property of molec¬ 
ular clouds, but instead just an artifact of limitations of the observational 
methods used to produce cloud surveys. Kegel (1989) and Scalo (1990) 
elaborated on this issue. Also, a number of observational studies have found 
significantly different scaling exponents for both Larson’s relations, or no 
clear scaling at all. For example, Carr (1987) has found a ^0.25 
M ~ , implying p ~ R~^'^ for 45 clumps in the Cep OB 3 molecular 

cloud. In fact, Carr notes that the clumps are not in virial equilibrium, and 
concludes that they must be expanding. Loren (1989) studied 89 clumps in 
the p Oph molecular cloud, finding p ~ R~^'^ and no a-R correlation. He 
points out that many of the clumps are extremely filamentary and that they 
may probably be the outcome of the passage of a shock wave. Falgarone et 
al. (1992) focused on low brightness areas of several molecular clouds, includ¬ 
ing both gravitationally bound and unbound structures (structures for which 
the actual masses are much smaller than the value obtained assuming virial 
equilibrium), finding a ~ R^''^, and a density-size plot with large scatter and 
an upper bound given by a Larson-type relation p ~ R~^. Vazquez-Semadeni 
et al. (1997) have found the same effect in numerical simulations of turbu¬ 
lent cloud formation in the ISM, in which in fact there is no unique p-R 
correlation. Plume et al. (1997), in a sample of maser-selected, massive star 
formation regions, find no significant a-R or p-R correlations either. 

Finally, Myers et al. (1995) have found an example of a cloud in which 
the kinetic and magnetic energies are comparable, yet the gravitational en¬ 
ergy is roughly 500 times smaller, indicating that self-gravity is negligible 
for this cloud. 


4. Discussion and conclusions 


From the discussion in the previous sections, a number of points can be 
concluded. 


1. “Virial equilibrium”, i.e., d?I/dt^ = 0 allows up to a constant rate of 
change in the moment of inertia I, even though it is commonly taken to 
mean a stationary equilibrium. 


2. Larson’s (1981) relations, eqs. (O) and (M), have been interpreted as 
a consequence of virial equilibriumgbetween self-gravity on the one hand 
and turbulent velocity dispersion and/or magnetic support on the other. 


^Mouschovias & Psaltis (1995) distinguish between magnetic virial equilibrium and mag¬ 
netic support, but the resulting balance equation is the same. 
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From the observational data discussed above, this interpretation may apply 
to strongly condensed clumps which, due to self-gravity, have decoupled 
from their surrounding medium, and thus require internal support against 
collapse. 

3. The above interpretation, however, is insufficient for explaining both scal¬ 
ing relations. One is still independent and requires an additional explanation. 
A critical mass-to-flux ratio has been invoked as the origin of the density-size 
relation assuming roughly constant magnetic field strengths (Shu et al. 1987; 
Myers & Goodman 1988b; Mouschovias & Psaltis 1995), but this assump¬ 
tion is not verified in the observational data. For example, the magnetic 
field strenghts discussed by Myers & Goodman (1988a) and Mouschovias 
and Pslatis (1995) span three orders of magnitude. Besides, examples of 
virialized clouds which however follow different scaling relations exist (e.g., 
Loren 1989; Fuller & Myers 1992; Gaselli & Myers 1995; Goodman et al. 
1997; see also Vazquez-Semadeni &: Gazol 1995 for a theoretical discussion). 
Planar shocks have also been invoked as an explanation of the apparently 
constant column densities (Larson 1981; Scalo 1987). 

4. It is possible that the density-size relation is not a true property of clouds, 
but only an observational effect due to the limited integration times used in 
surveys (Larson 1981; Kegel 1989; Scalo 1990), which tend to select constant- 
column density obejcts. In this case, surveys using larger integration times 
would exhibit an increasingly larger scatter in density-size plots, as in the 
data of Falgarone et al. (1992) or the numerical data of Vazquez-Semadeni 
et al. (1997), the latter being free from such limitations. In fact, large scat¬ 
ter was already present in the Myers & Goodman (1988b) data. Although 
they interpreted it as uncertainty in the data rather than real scatter about 
virialization, the latter possibility is equally feasible. A study that seemed to 
hnd exceptionally constant column densities while claiming a large dynamic 
range (Wood et al. 1994) has been questioned by Vazquez-Semadeni et al. 
(1997). 

5. The velocity dispersion-size relation does not appear subject to the pos¬ 
sibility of a spurious origin, and therefore it is probably real when detected, 
although in many cases it is not (e.g., Loren 1989; Plume et al. 1997). For 
the cases when it is present, a number of physical mechanisms are plausible 
candidates. If the density-size relation is not verified, then the standard ar¬ 
guments based on virial equilibrium between gravity and velocity dispersion 
cannot be invoked (recall the a-R relation is a consequence of virial equi¬ 
librium and the density-size relation). In particular, the possibility that the 
a-R relation is a consequence of the characteristic energy spectrum of an 
ensemble of shocks in the flow has been suggested by a number of authors 
(e.g., Passot et al. 1988; Padoan 1995; Gammie & Ostriker 1996; Fleck 1996; 
Vazquez-Semadeni et al. 1997). In this case, the origin of the dispersion-size 
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relation would be completely independent of the virial condition, explaining 
why the relation can be present even in the absence of a density-size relation. 

6. In summary, the Larson’s relations and their virial equilibrium interpre¬ 
tation probably apply to relaxed, strongly self-gravitating clumps. Regions 
which do not exhibit clear scaling relations often include either massive cores 
(Caselli &: Myers 1995; Plume et ah 1997) or regions with strong evidence 
of recent perturbations (Loren 1989). In fact, the dense cores studied by 
Plume et ah were selected by the presence of H 2 O masers, suggesting strong 
excitation mechanisms as well. Such “perturbed” regions are likely not in 
virial equilibrium, thus also being transients (e.g., Magnani et ah 1993), or 
at least having strongly fluctuating moments of inertia (shapes and internal 
mass distributions). 

Acknowledgments. The author wishes to thank A. Goodman, T. Pas- 
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